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O ■ Abstract 

o 

f~^ I We present model predictions for the spectrum of CP in a periodic box and 

C^ ' use them to interpret the strong finite size effects observed in lattice simulations at 

medium values of N. The asymptotic scaling behaviour of alternative lattice actions 
j^ I is discussed along with some aspects of multigrid algorithm efficiency. 



Introduction There has been considerable recent interest in simulations of the lattice 
(jpN-i model [jl], 0, ^, |, §. Advances in the development of non-local Monte Carlo algo- 
rithms have given added impetus to studies of the model at large A^ where the spectrum 
remains unconfirmed by analytic or numerical analysis. Conflicting evidence has been pre- 
sented for perturbative (two-loop) scaling at currently accessible values of the bare coupling. 

In this letter we analyse some key factors involved in such studies: the dependence on 
lattice action, the efficiency of the algorithm and, most crucially, contamination of finite size 
effects which are particularly significant in the present model at medium and large A^. We 
demonstrate the latter effect by means of a simplified analytic model and present some new 
data. 

The simplest lattice action for the model containing an explicit gauge field is 

S = P/2Y: [zKx)Uix,y)z,iy) + h.c] (1) 

<xy> 

where (3 = 2N/g'^ and z(x) is a complex A^-vector with unit norm. U{x, y) is a gauge 
connection between nearest neighbour sites x,y and is an element of U{1). This action, 
which was used in the studies P, ^, corresponds in the classical continuum limit to the 
Lagrangian 

^ = ^(^.z)t . {D,z) (2) 

where the gauge field appearing in the covariant derivative may be related to the 'matter' 
field z{x) via its equation of motion. For computational convenience, the lattice gauge 
variable can be integrated out to obtain an effective 2D spin model 



whose action is 



The alternative action 



Z = f[dz][dz*] exp{-(3S{z)} (3) 

I3S = -Y: \nIo{f3\z*{x)-z{y)\) (4) 



<xy> 



Sq = -2kY: \^*{x)-z{y)\' (5) 

<xy> 

can be shown to have the same classical continuum limit as (|l]) with the identification of 
bare couplings k = N/2g^. This latter action was studied by |]I|, ^, Q. 

At large A^, an effective Lagrangian may be derived [^ 0] which corresponds to A^ identical 
charged scalar fields subject to a (confining) ^(1) gauge potential 

Viv)=gj\r\ (6) 

where the effective coupling is related to the mass of the charged scalar rrif 



9f = ^/y ]^ (7) 

In the large A^ limit, the low-lying spectrum can be estimated in a non-relativistic approx- 
imation via a Schrodinger equation for the two-body relative motion in the potential V(r). 
1^. For medium A^ (4 ... 8), a previous lattice analysis showed that the predictions of this 
simple model are at best unreliable 0. The lowest state observed was indeed the adjoint, 
but there was little sign of the next predicted state, an SU{N) singlet. It was concluded 
that either the lattice singlet operators used were inefficient or that the breakdown of the 
non-relativistic model approximation at these A^ values was to blame. 



Model for finite size effects All analyses of CP^~^ have noted strong finite-size effects. 
Detailed knowledge of the A^-dependent potential allows one to estimate the size of lattice 
required to be free of such effects 0. Most analyses have required L/^ > 6 for this reason. 
In fact, one can detect the effects of working in a periodic box even above this limit especially 
at large values of N where the flat potential leads to a weakly confined sytem and hence 
broad wave functions 0, |^. 

It is interesting to use the (albeit imperfect) non-relativistic model described in the first 
section to investigate these effects more quantitatively. The potential implied by eqns (^, |^) 
when made periodic with period L gives rise to a band rather than discrete spectrum (see 



any Quantum Mechanics text for the treatment of the Kronig-Penny potential |Tll). The 
discrete eigenvalues determined by the zeroes of Airy functions [§ give way to allowed bands 
whose limits are determined by combinations of Airy functions 

AA cos{kL) = [B'A + A'B) {B'A + A'B) - 2AA'BB' - 2BB'AA' . (8) 

Here, k is some real parameter and L is the period of the potential: the spatial extent of the 
lattice in our application. Also 

A = B'A- A'B, A = B'A- A'B (9) 

where A = Ai{—b) and B = Bi{—b) are the standard Airy functions, with A', B' the 
corresponding derivatives. Similarly A = Ai{aL/2 — b), and so on, where the parameters a, 
b are related to the potential slope 2Vo/L, the required band energy E and the reduced mass 
of the two-particle system /i = mf/2: 

. = (M£)V3, 6/a=^ . (10) 

In the present application, Vq = g'jL/2. All scales in the model are set by the confined, and 
so unobservable, fundamental mass ruf. To compare predictions with data (see below) the 
relationship to the physical bound state mass must be used. 

In Fig 1, is shown the band spectra for A^ = 4 and 8 in units of rrif as a function of 1/L. 
These have been checked against the results of numerical solutions of the corresponding 
Sturm-Liouville problem. Using the numerical approach we have also been able to replace 
the model potential (^ by that measured in lattice simulations [0]. The qualitative results 
remained similar to those seen in Fig 1. The spectrum at 1/L = agrees with that of Haber 
et al. 1^. i.e. a low-lying positive parity adjoint state followed by well separated negative 
parity and further excited states. As the size of the periodic box decreases in physical units, a 
band structure develops implying that zero momentum correlators in lattice simulations will 
no longer be dominated by a single isolated state. At increasing Euclidean time separation, 
the contribution from the lower levels within the lowest band will dominate so that the 
effective mass me//(t) defined by 

C(t) = const. X [exp(— me//t) + exp(— me//(-L — t))] (11) 

shows a drop-off in t rather than the plateau characteristic of an isolated state. One may 
construct simple models for the overlap of the band states in Fig 1 with the measurement 
operator employed in C{t) which demonstrate the transition from meff{t) ~ mo at moderate 
Euclidean time separations to rrieffit) ^ mo — 6/2 at large separations. Here, mg is the centre 
of the energy band for some lattice size and S its width and all units set by my = 1. 



Since the non-relativistic model is known to be quantitativley unreliable [Q it is difficult to 
make other than qualitative statements on the implications of this effect. As an example, 
for A^ = 4 and k = 2.9 (action (^), ^/a = 18.5. According to the non-relativistic treatment 
§ for AT = 4 

m = 1/^= (2 + 4.54)m/ (12) 

Thus, a lattice of size L/a = 128, say, corresponds to L = 1.06mJ^. According to Fig 1, 
this implies a very large broadening effect indeed. Of course, we know that the system is 
actually relativistic since eq . (|12D implies v/c ~ 0.95 so model predictions are unreliable. On 
the other hand, the string tension a implied by the large A^ model (cr/m^ = 37r/6.54^ = 0.22) 
is comparable with that measured (0.34(2)) [|2| when expressed as a dimensionless ratio. 

Hasenbusch and Meyer |^ have presented correlator data at very large time separation, 
which the authors claim is essential to obtain the true mass-gap in the presence of other 
states. While this is in general true, the above model demonstrates that for systems bound 
loosely by a shallow potential such as large A^ CP^~^ the effect of periodic boundaries 
could be such that large time separations isolate the wrong state i.e. not that which is 
physically relevant in the bulk, continuum limit. In the absence of detailed published data 
for nieffit) we have investigated this point by studying CP^ and CP^ at weak coupling 
on small (64^) and large (128^) lattices. Fig 2 shows typical results. There is some small 
but significant effect apparent, particularly for A^ = 8 where the potential is flatter. Note 
that the errors shown are statistically correct at each time separation but are correlated 
in Euclidean time. Clearly the effect is considerably smaller than the naive non-relativistic 
model (Fig 1) predicts. Nonetheless, we would suggest that there is a significant source of 
systematic error which must be allowed for in future high statistics analyses. To be free of 
finite size effects, it would seem prudent to study systems for which [Q 

L/e > cN (13) 

where c is of the order of 2. 

For much smaller lattices, one sees a reversal to the more usual finite size effect of raised 
effective masses (reduced correlation lengths) 0. 

Scaling properties Analyses using the quadratic action (^ have shown conflicting evi- 
dence for scaling of the mass gap (adjoint state) [|^, ^. The first study with relatively low 
statistics and covering the lattice correlation range ^/a < 30 showed some evidence for scal- 
ing (constant ratio of dimensionful quantities) and for asymptotic scaling of the mass gap 
itself I, [11 

ma = (m/Ajv)(^)^ exp{ -} (14) 

r 9 

This analysis [Q] relied largely on a local Metropolis algorithm. However, a multigrid ap- 
proach enabled Hasenbusch and Meyer [|] to achieve larger statistics and to probe larger 
correlation lengths (< 77). Clear evidence of non-asymptotic scaling was presented. 

Irving and Michael used the gauge-explicit action (|l|,|D for the models with 4 < A^ < 8 and, 
using a hybrid local/cluster algorithm, presented data showing evidence of scaling (ratio of 
topological charge and mass-gap) and for asymptotic scaling of the mass-gap for g^ < 1.1. 
A subsequent analysis [^] using a hybrid over-relaxation/heat bath algorithm on the same 
action with explicit gauge fields has provided further convincing evidence that both scaling 



Quadratic action 


Explicit gauge action 


1/^^ 


ila 


mjk 


1/^' 


e/« 


m/A 


1.25 


4.48(8) 


514(9) 


0.909 


3.1(1) 


102(3) 


1.35 


8.7(3) 


478(12) 


0.977 


4.9(1) 


97(2) 


1.4 


12.8(3) 


436(11) 


1.000 


5.3(1) 


98(3) 


1.45 


18.5(3) 


406(8) 


1.111 


10.8(3) 


95(3) 


1.5 


30.0(2) 


337(15) 


1.190 


18.5(7) 


88(3) 



Table 1: The scaled adjoint mass m/A in lattice CP^ using the quadratic action and the 
explicit gauge action (see text). Also shown are the corresponding hare couplings and corre- 
lation lengths in lattice units (^ = 1/m). 



and asymptotic scaling properties are superior for this action when compared to the quadratic 
action (0). 

We have used a variant of the multigrid algorithm of Hasenbusch and Meyer [^ to compare 
the asjTiiptotic scaling in the A^ = 4 model (|14|) for both actions using the same analysis tech- 
niques. The t-dependent effective mass-gap was extracted from zero-momentum correlators 
The errors on the t-independent values (see discussion above) were extracted 



11 



usmg eqn. 

using bootstrap techniques. Typically 20000 multigrid V-cycles with measurements every 
second cycle were made at each coupling. As can be seen from table |I], the corresponding 
values of m/A are not constant in 1/g^ but are significantly more slowly varying for the 
explicit gauge action S. If one parameterises the the scaling violation by 

m/A ~ i^/a)-'^ (15) 

then the data presented in table |l| yield Zg = .18 ± .02 (.07 ± .02) for the quadratic (explicit 
gauge) action. 



Algorithms The data presented above was obtained using an adaptation of the multigrid 
algorithm proposed by Hasenbusch and Meyer ||^. The particular improvement tested 
was to use random SU{2) subgroups rather than f/(l) and to employ a different subgroup 
and rotation at each level of the V-cycle. Using this, we made careful measurements of 
the decorrelation dynamical exponent z using the spin susceptibility as a good measure for 
the self-consistent autocorelation time. For example, we found values close to 1.0 for the 
quadratic action (Kl) 



Z{N=4) 



1.05 ±.08, 



Z{N=8) 



1.12 ±.10 



(16) 



i.e. considerably larger than those reported by ^ but comparable with those reported 
by Wolff 1^ who used an overrelaxation algorithm. When we used a fixed U{1) rotation 
for each block within the V-cycle, in the spirit of |Q the decorrelation performance of 
our algorithm was degraded. The multigrid performance on the action (^ was similar: the 
dynamical exponent for CP^ was 1.00 ±.07. Of course, all these algorithms offer considerable 
improvement over the local Metropolis, single cluster, and hybrid algorithms used in earlier 
work MM. 



Conclusions We have presented a model for non-trivial finite size effects in CP^^^ and 
new data relevant to this, to the question of precocious scaling and to the efficiency of 
multigrid algorithms. Effects ascribable to periodicity as predicted by the simple model 
have been demonstrated but at a much lower level. We have confirmed earlier evidence 
10, ^ that the explicit gauge action for lattice CP^~^ has superior scaling properties to the 
more commonly used quadratic action. Finally, we reported measurements of the dynamical 
exponent for a multigrid algorithm which represent a considerable improvement on local 
algorithms but which fall short of other reported data for similar algorithms [^ . 
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Figure captions 

1. The low- lying part of the predicted energy spectrum due to a periodic potential. 
The energy bands allowed by (^ are indicated by the shaded regions bounded by 
solid(dashed) lines for A^ = 4(8). Units are set by wie// = 1. 

2. The effective mass defined by eqn. |TT] for (a) A^ = 4 at k = 2.9 and (b) A^ = 8 at 
K = 5.5 The crosses (circles) denote data on 128^(64^) lattices. The solid line is the 
asymptotic mass gap used in the scaling analyses. 



